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^— ' The mean first passage time (MFPT) for pliotoexcitations diffusion in a funneiing potentiai of 

, artificiai tree-like liglit-liarvesting antennae (pfienylacetylene dendrimers witli generation-dependent 

segment lengtlis) is computed. Effects of tfie non-linearity of the realistic funneling potential 
and slow random solvent fluctuations considerably slow down the center-bound diffusion beyond 
^ , a temperature-dependent optimal size. Diffusion on a disordered Cayley tree with a linear poten- 

^ ' tial is investigated analytically. At low temperatures we predict a phase in which the MFPT is 

dominated by a few paths. 

00 



O 
(N 



s 

I 



o : 
o 



Dendrimers constitute a new class of nanomaterials with unusual tree-like geometry and interesting chemical, 
transport, and optical properties |^-|^. Two families of Phenylacetylene dendrimers have received considerable recent 



^ , attention. In the compact family, the length of the linear segments is fixed, whereas in the extended family it 



increases towards the center creating an energy funnel in that direction (Fig. 1). This latter family may, therefore, 
serve as artificial light-harvesting antennas, as has been demonstrated experimentally It has been conjectured by 
Kopelman et. al., based on optical absorption spectra [0, that electronic excitations in these dendrimers are localized 
on the linear segments. This has been confirmed in theoretical studies which showed that the relative motion of 
photogenerated electron-hole pairs is confined to the various segments and energy-transfer may then be described by 
'"^ , the Frenkel exciton model. The time it takes for an excitation that starts at the periphery to reach the center, and 
^ ■ its dependence on the molecular size (number of generations g) and the funneling force, were calculated. The latter 
O results from the interplay of entropic (or geometric, i.e. the branching ratio c = 2) and energetic factors 
I ^ I These pioneering studies, however, assumed the funneling force to be constant since the energy e{n) varied linearly 

with n. {n=l,2,. . . ,g in the generation number, at which the segment length is ^ = — rt + 1 monomers. See Fig. 
' ' 1.) In addition, e{n) were assumed to be fixed. In reality, interactions with other degrees of freedom (solvent and 
intramolecular vibrations), induce fluctuations in e(n) which may span many different timescales. Here we consider 
slow (quenched) fluctuations compared with the exciton trapping times which are typically in the picosecond range 
IHjH- Nonlinear j|] and single molecule [|lo| spectroscopy in liquids, glasses and proteins typically show nanosecond to 
millisecond bath motions responsible for spectral diffusion. Slow vibrational motions that can be treated as static dis- 
\ order dominate the photoinduced energy transfer dynamics of photosynthetic antenna complexes . Fast (annealed) 
fluctuations do not change the behavior qualitatively. In addition, we explore different degrees of correlations among 
the various energy ffuctuations. In the absence of correlations, we obtain the standard diagonal disorder (random 
' energy) model. If the energy differences of nearby segments are coupled to independent baths, we obtain a random 
^ ' force model. Other types of correlations are possible as well. In a recent paper Bar-Haim and Klafter extended 
I their studies to a special kind of disorder with unique correlations between the hopping rates. They also looked at 
' O ■ the effect of including a single impurity. 

g ■ In the present Letter, we address two crucial properties of the funneling potential of dendrimers: (i) The nonlinear 
^ ' dependence of the exciton energy on n is taken from the electronic structure calculations of Tretiak, Chernyak, and 
1^ Mukamel iQ. (ii) The effects of general and realistic types of quenched disorder are investigated. Our study has direct 
. ^ implications on the design of the dendrimers with improved light-harvesting efficiency. We find that the nonlinear 
potential drastically reduces this efficiency for dendrimers larger than a polymer-specific optimal size (which is found 
H I to be in the range of 7-9 generations at room temperature). At the same time the effect of disorder in slowing down 
■ - - ■ the excitonic migration towards the trap, is much more severe for dendrimers larger than this optimal size. Thus, on 
both counts, efforts to make the dendrimers larger than this size will provide only a modest pay-off in the number of 
excitons reaching the active center (even though the total photon adsorbance will increase [Q). 

Our results differ drastically from the ones found for a linear potential. However, they can be qualitatively under- 
stood in the context of the different regimes of diffusion in such a potential. Therefore, we first study the effect of 
disorder on the diffusion in a linear potential. Interestingly, we find a dynamic transition into a highly disordered 
phase akin to the equilibrium replica symmetry breaking found in other random statistical models [|l3| . 

The theoretical investigations focused on the mean first passage time (MFPT) , which provides an adequate measure 
for the efficiency of the trapping. Another quantity of interest is the mean residence time (MRT) , the average time 
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spent on a site of the tree (i.e., on a segment of the dendrimer). The starting point for the analysis are coupled master 
equations for the probability p{i, t) of the exciton to reside at time t on the «th segment. If the potential only depends 
on the generation n, the problem reduces to an effectively one-dimensional one for p(n, t), which obeys the equations: 
Pn{t) = Tnj^ipn+iit) + Rn-iPn-i{t) — {Tn + Rn)Pn{t) , whcrc T (i?) are the rate constants for transfer towards the 
periphery (center) and absorbing (reflecting) boundary condition are used at n = {n = g). 

We define the detailed-balance ratio ^(n) — i?„/r„+i = cexp{— /3[e(n+l)) — e(n)]}. For simphcity we assume r„ = 1 
for all n. For excitation starting at the periphery, the MRT on the nth site is given by 0: tn — {J2^m=i TljCn-m + 
The MFPT is given by the sum over the tnS: Tg = X]n=i^"- ^'-^^ ^ linear potential ^(n) = = cexp(— /?/), where 
/ is the n-independent potential-difference and P is the inverse temperature (room temperature is assigned unless 
specified otherwise), one finds 0: 

r,=eO(^^-^, Co^l. (1) 

Eq. (1) shows three regimes for the scaling of Tg with g: For > 1: Tg ^ exp{glTi^Q) (exponential regime) for i^p < 1^ 
Tg ^ g (linear regime) and at the transition, = 1 (for a critical force fc = Inc/P), the behavior is purely diffusive 
(quadratic): Tg ^ g'^ 0. 

Nonlinear TDHF potential. 



The exciton energies of linear segments of acetylene units have been calculated in |12| using the time-dependent 



Hartree Fock (TDHF) technique. These energies can be fitted to the following nonlinear expression |14 



e(n) = A 



n 0.5 



g-{n- 1) 



(2) 



with A=2.80 ±0.02 (eV) and L=0.669 ±0.034. 

The MFPT for an exciton generated at the periphery on this potential is depicted in Fig. 2, for different tempera- 
tures. Its variation with g differs substantially from that of a linear potential. Although the MFPT depends linearly 
on g for the first few generations, it gradually crosses over (for n > 7) to an exponential behavior with increasing g. 
As the temperature is reduced, this crossover scale increases due to the entropic effect. 

The MRT is not a monotonic function of n. For large g, at the generations near the periphery (g — n <^ g) the 
tunneling force is strong and overcomes the entropic effect of c > 2 which "pushes" the exciton outward. Near the 
center, however, the larger is g, the weaker is the tunneling force and the entropic term dominates. As a result the 
exciton spends most of its time at some intermediate generation n = n* (g) close to where the two effects nearly cancel 
and the net thermodynamic force is minimal. In Fig. (3), we display the reduced free-energy u(n) = (3c{n) — nine, 
and the MRT, vs n for g ^ 14. The MRT is maximal at n*(14) = 8 close to where u{n) is minimal. As g increases 
beyond g ~ 9, both n* (Fig. (3)) and the energy difference Au{n*) = u{l) — u{n*) increase with g. Hence, the time 
to reach the center from n* grows exponentially with g (the time to arrive at n* is much shorter). 

Random Cayley tree with a linear potential. 

We begin the investigation of the effects of disorder by first considering a linear potential. Four models were studied: 
(i) Random intergenerational energy - The same random energy is assigned to all the segments in a given generation. 
(a) Random intergenerational force - The energy differences between consecutive generations are randomly distributed. 
(Hi) Random intersegment energy - All segment energies are random and uncorrelated, and (iv) Random intersegment 
force - All energy differences between neighboring segments are random. The mapping to an effective one-dimensional 
model applies for the first two models only. We set ^„ = exp(— /3Ae„) with Ae„ — e„+i — e„. The behavior for each 
of the models is as follows: 

(i) Random intergenerational energy: This model whereby e(n) have an identical distribution P(e) (with < e„ >= 0, 
<> denotes average over the disorder), has been investigated ||l5| with the following implications to our system: 
We define 77 = exp{—(3e) and assume finite values for < 77^^ > and < 77^^ >. The average MFPT (up to g- 
independent constants) is given by : < T{g) >=< 77 >< rj^^ > TQ^g) with the three regimes of behavior (linear, 
quadratic and exponential) depending on ^0 being <,=,or > 1. For example, for a Gaussian distribution Pcif-) = 
(l/\/27rA) exp{— e^/2A^}, < T{g) > is enhanced by a ea;p(/3^A^) factor compared with To{g)- The disorder then 
slows down the tunneling, but leaves the g-dependence unchanged. Another important effect of disorder is to induce 
fluctuations in the MFPT < Ar^ >=< (r- < r >. The relative fluctuations 5t/t = < Ar^ >i/2/< r > scale as 

Wff' linear regime and at the quadratic point, and is a ^-independent constant in the exponential regime. 

(ii) Random intergenerational force: This model assumes Ae„ to be independently distributed according to a 
distribution P(Ae) and is reviewed in Un. Proceeding along similar lines, we find the following regimes: (a) < ^ >< 1: 
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< T{g) > is linear in g. < T'^{g) >~ g^ as long as < >< 1, while if < >> 1, < r'^ig) > grows exponentially 
with g as the typical behavior begins to differ from the average. (6) < ^ >> 1, but < In^ >< 0: < r((?) > is 
exponential in g and so are the relative fluctuation 6t/t. The average is determined by rare configurations, while 
the more representative Ttyp = exp < \ogT{g) > behaves as g" with a — < S{ln^)'^ >/2 < In^ > [— /3^A^/2(ln^o) for 
PG(Ae)]. (c) < In^ >= |l^: < T{g) > is exponential in g while Ttyp{g) is exponential in ^ {d) < In^ >> 
both < T{g) > and Ttyp{g) diverge exponentially with g. 

(Hi) Random intersegment energy: We extended the calculations of model (i) to the fully random-energy tree 
with almost identical results. The main difference is that, for intersegment disorder, additional fluctuations in the 
MFPT arise from distinct initial sites at the periphery. In all regimes their g dependence is identical to that due to 
denderimer-to-dendrimer fluctuations. St/t (from both effects) saturates, at large disorder, to a value smaller than 
the corresponding one in model (i). 

(iv) Random intersegment force: The techniques applied in model (ii) may not be simply generalized to the fully 
random-force tree. Using the replica trick, we have found a one-step replica symmetry breaking (IRSB) transition 
[^3| between a weakly disordered (or high-temperature) phase in which all paths contribute to T{g) and a highly- 
disordered (low-temperature) phase in which only a small number of them do. Mathematically, the transition is 
determined by the parameter < m < 1 for which the value of [^o < exp(— /3mAen) >]^/'" is minimal, m = 1 in 
the high-temperature phase, limited to /3A < y/2\nc for Po{Ae), for which /c = lnc//3 -|- /3A^/2 in this phase. In the 
low-temperature "glassy" phase m — /3A/\/21nc < 1, and the critical force fc = A-\/2 Inc is temperature-independent. 

Nonlinear TDHF potential with quenched random energies. 

We have carried Monte Carlo simulations of the effect of disorder on the MFPT for the TDHF nonlinear potential. 
Disorder is introduced by allowing the energies e(n) to fluctuate uniformly in the range of ±2% around their pure 
TDHF value eo(n) (Eq. (2)) 

Fig. 4 shows the disorder-averaged (over 10* realizations) MFPT for both intergenerational and intersegment types 
of disorder [such as models {i) and {Hi)], compared with the pure TDHF potential. The effect of disorder to increase 

< T{g) >, is clearly more pronoimced in the exponential regime. The plots of < r(g) > for intergenerational and 
intersegment disorder are indistinguishable. Their relative fluctuations, however, do differ for large g. 6t/t saturates 
to two different values (inset of Fig. 4) that of the intergeneration disorder being larger. All these results may be 
understood from our analytical analysis of models (i) and {Hi) above. 

In summary, we have shown that the realistic non-linear potential induces an effective tunneling only for molecules 
smaller than some optimal size. For larger dendrimers the free energy has its minima at n*{g) where the excitons 
spend the largest amount of time. For a large production rate of long-lived photoexcitations, we expect the excitons 
to accumulate at n*{g), at which case the single-exciton picture is not applicable and the exciton-exciton interactions 
and annihilation processes need to be accounted for Disorder slows the excitation diffusion towards the center. 
This effect is always stronger if the MFPT of the corresponding pure system is in the exponential regime. Hence, in 
order to minimize this slowing down due to randomness, the dendrimers have to be smaller than the same optimal size. 
Finally, we have uncovered a new effect of increasing force-randomness (or lowering temperature) on the dynamics: 
the MFPT is dominated by a few paths along the tree. It should be pointed out that not only the average A4FPT but 
the complete distribution of T{g) is readily available experimentally from the time resolved fluorescence profile of the 
antenna or an acceptor at the center. It will be interesting to explore experimentally this distribution and especially 
its dependence on the number of generations, and the temperature. 
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FIG. 1. Extended Phenyacetylene Dendrimer with g = 4: generations. 



FIG. 2. The MFPT vs g for the TDHF nonlinear potential at different temperatures (Inset: same for g < 5 magnified). 



FIG. 3. Left: The MRT vs n for the TDHF potential, for different values of g. Right: The reduced free-energy u{n) for 
5 = 14. 

FIG. 4. The MFPT for the TDHF potential with both types of random energy (indiscernible from each other) compared to 
that of the pure system {Inset: Their relative rms variations vs g). 
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